Quantum interference in macroscopic crystals of no n- metallic Bi2Se3 
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Photoemission experiments have shown that Bi2Se3 is a topological insulator. By controlled 
doping, we have obtained crystals of Bi2Se3 with non-metallic conduction. At low temperatures, we 
uncover a novel type of magnetofingerprint signal which involves the spin degrees of freedom. Given 
the mm-sized crystals, the observed amplitude is 200-500 x larger than expected from universal 
conductance fluctuations. The results point to very long phase breaking lengths in an unusual 
conductance channel in these non-metallic samples. We discuss the nature of the in-gap conducting 
states and their relation to the topological surface states. 

PACS numbers: 72.15.Rn,73.25.+i,71.70.Ej,03.65.Vf 



A new class of insulators with non-trivial topological 
surface states has been predicted P, 0, 0, 0, Q ■ The sur- 
face states of these topological insulators are chiral, and 
protected from disorder by a large spin-orbit interaction 
that aligns the spin transverse to the wavevector. Angle- 
resolved photoemission spectroscopy (ARPES) has been 
used to detect these surface states (SS) in Bii_a;Sba; [1]. 
The spin polarization of the SS in Sb was confirmed by 
spm-r esolved ARPES [1]. More recently, ARPES experi- 
ments showed that Bi2Se3 [lOiliii and Bi2Te3 [l^ are the 
simplest topological insulators, with only a single Dirac 
surface state. 

Due to its large band gap (300 mV), Bi2Se3 is a very 
attractive platform for exploring the transport proper- 
ties of the topological states. However, as-grown crystals 
of Bi2Se3 invariably display a metallic resistivity profile, 
with the bulk chemical potential fii, pinned to the con- 
duction band. By chemical doping, we have obtained 
crystals of Bi2Se3 in which fii, falls inside the gap. Al- 
though the residual conductance at low temperature T 
(0.3 K) appears to be still dominated by bulk conduc- 
tion channels, we observe a novel type of conductance 
fluctuation phenomenon. In a swept magnetic field H, 
the fluctuations retrace reproducibly, analogous to the 
magnetofingerprint signal reminiscent of universal con- 
ductance fluctuations (UCF) [l^. However, the finger- 
print signal here is observed in mm-sized bulk crystals 
rather than mesoscopic samples. Moreover, field-tilt ex- 
periments show that the carrier spin plays a role in gen- 
erating the flucutations. We discuss the highly unusual 
nature of the large magnetofingerprint signal and possible 
connections to lattice dislocations that bear topologically 
protected states. 

The ARPES results reveal that Bi2Se3 has a single 
Dirac surface state (SS) that crosses the bulk energy gap 
(Fig. [1^). In as-grown crystals, electrons donated by Se 
vacancies pin /i^ to the conduction-band (CB) edge. Re- 
cently, Hor et al. showed that doping with Ca converts 
the crystals to p-type conductors [14[ . By tuning the Ca 
content x in Ca2,Bi2_a;Se3, we progressively shift from 
the CB to inside the gap, and then into the valence band 
(VB). Samples with /i;, in the CB or VB display Shub- 
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FIG. 1: (a) Resistivity p vs. T in 4 samples (G3-G6) of 
Ca2,Bi2-iSe3 lightly doped with Ca to bring Hb into the gap. 
Samples with ^i, not inside the gap (M2, M5 and MIO) dis- 
play a metallic T dependence (shown x20). The inset is a 
sketch of the surface states [T^l crossing the gap from the VB 
to the CB. (b) The LL index plot vs. field minima in the 
SdH oscillations observed in 8 metallic samples (Ml- ■ • Mil). 
Negative (positive) index n represents the electron (hole) FS 
pocket. As Hb is lowered from the CB to VB, the FS area 
Sp ~l/slope decreases before rising again. 



nikov de Haas (SdH) oscillations which we have used to 
measure the caliper area Sp of the bulk Fermi surface 
(FS). The carrier sign was found by the thcrmopower 
and Hall effect. The index plots of the Landau Levels 
(Fig. [TId) shows that Sp decreases (Samples M1^M3) 
as enters the gap from the CB, and then increases as 
jib exits the gap and moves into the VB. The metallic 
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FIG. 2: Curves of p vs. H in Sample G4 at 0.3< T <8 K. 
The curves are plotted vs. H (Panel a) and vs. InH (b). The 
MR displays a sharp anomaly of amplitude 10%-R in weak H . 
In addition, large conductance fluctuations of amplitude up 
to 0.5% are resolved. In Panel (b), the plot of p vs. log_ff at 
0.3 K shows that p{H) ^ \nH over 2 decades in H (dashed 
line). The slope is equivalent to dG/dlnH = 200e^//i. The 
correlation function C{B,T) of the fingerprint signals (H||z) 
is plotted in Panel (c) for T = 0.3 K (bold curve), 0.85 K 
(medium) and 1.8 K (thin). C{B) oscillates vs. H instead of 
decaying as a power law. 



profiles of p are shown in expanded scale in Fig. [T^ (M2, 
M5 and MIO). Using the changes in Sp to guide the dop- 
ing, we have obtained non-metallic crystals in a narrow 
doping window 0.002< x <0.0025, with p,i, lying inside 
the energy gap. SdH oscillations are not resolved in any 
of the non-metallic samples. 

As plotted in Fig. [T^, p of the non-metallic crystals 
(G3-G8) undergoes an increase to very large values (20- 
100 mrJcm) as T falls below 130 K. Below 20 K, p ap- 
proaches saturation instead of diverging as in a semicon- 
ductor. Although the conductances G of the non-metallic 
crystals are very poor (Table|T|, they are still 1000-8000 x 
the universal conductance e'^/h (e is the charge and h is 
Planck's constant). As we will describe, draining away 
the high-mobility electrons in the CB uncovers a conduc- 
tance channel of a highly unusual kind. 

The magnetoresistance (MR) in Sample G4, measured 
with H||c||z (with the current I||x), is displayed in Fig. 
[5^ for T ~ 0.3-8 K. In each curve, the most promi- 
nent feature is the pronounced weak-field anomaly, which 
deepens to a sharp cusp at iJ = 0. Near the cusp, p vs. H 
follows a logarithmic behavior extending over 2 decades 
in H at 0.3 K (Fig. HJa). Moreover, at low T, large 
conductance fluctuations (~ 0.5% p) are apparent. The 
retraceability of the fluctuations versus H distinguishes 
them from random noise (the correlation function shown 
in Panel c is discussed later). 

Traces of the conductance fluctuations SG{T, H) — 
G{T, H)—Go{T, H), relative to the smoothed background 
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FIG. 3: Curves of the magneto-fingerprint signal SG{T, H) vs 
H in Sample G4. The field H is aligned with z (in Panel a), 
with y (in b), and with x||I (in c). In each panel, curves are 
shown for five T between 0.3 and 4 K (in ascending order). 
For clarity, adjacent curves are displaced vertically by 10 e^//i. 
At 0.3 K, both up-sweep and down-sweep traces are shown 
superposed to emphasize retraceability. 



Go (T, H) , are displayed in Fig. [3] for 3 field orientations 
(H||z, H||y and H||x in Panels a, b and c, respectively). 
In each panel, we have emphasized the retraceability of 
SG{T, H) by superposing the up- and down-sweep curves 
at 0.3 K. Significantly, they remain large even with H||I 
(Fig. [5J:). We have followed the fluctuations to fields of 
32 T in G4 and G8. The root-mean-square (rms) ampli- 
tude is nearly unchanged between 4 and 32 T (by con- 
trast, SdH amplitudes should grow exponentially). In all 
non-metallic crystals studied at T <3 K (Table [!}, the 
magnetofingerprint is present. 

The existence of such large conductance fluctuations, 
with amplitudes ±10 e^/h, is remarkable in a bulk 
mm-sized crystal (2x2x0.05 mm^ for G4). Because 
magnetofingerprints imply field modulation of the inter- 
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TABLE I: Sample parameters, c is the crystal thickness along 
c. Values of G and rms5G (in /h) and p are measured at 0.3 
K (except for G3, which was not cooled below 4 K). We define 
A = Aorb + Aspin (scc Eq. [2]). The crystals are of nominal 
size 2 mmx 2 mm x c. The Hall density uh ~ l/eRn is 
inferred from the Hall coefficient Rh- 
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ference between conductance contributions from many 
channels, it is natural to compare them with Aharanov 
Bohm (AB) oscillations and UCF investigated in meso- 
scopic samples. However, we will argue that they belong 
in a new category. 

To see this, we recall the main features of the AB os- 
cillations and UCF. In mesoscopic, multiply-connected 
samples (single loops or arrays of loops), two types of 
conductance oscillations are observed with periods in flux 
(j) = h/2e and h/e, respectively. The former, called AAS 
oscillations [TgI . \V1\ , involves interference between a pair 
of time-reversed states that circumscribe a loop in op- 
posite directions, and are observed only in very weak H 
(<100 Oe). The latter (AB oscillations), arising from in- 
terference between waves that traverse opposite arms of 
the loop, survive to very large H (>10 T). However, the 
AB oscillations do not ensemble average. When the sam- 
ple size L exceeds the phase-breaking length ~ 1 /im, 
the amplitudes decrease as ^ L^jL jl9j. By contrast, the 
AAS oscillations have been observed - with weak if - in 
giant arrays (10^ loops) mm in size [T8| . 

Apart from periodic oscillations, there exist weak, ape- 
riodic fluctuations of the conductance G vs. H (UCF) in 
simply-connected mesoscopic samples [1^ US] . For elec- 
trons undergoing quantum diffusion in a phase-coherent 
region (L ~ L^), the fluctuations dG do not self-average. 
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The rms value rms[(5G] is ~ eV/i [13|, |20|, |2l|, I23 . Mod 
ulation of the interference by H results in the mag- 
netofingerprint trace. The amplitude of UCF is also 
sharply suppressed if L exceeds L,^. In metals, the phase- 
coherent region is actually cut off by the thermal length 
Lt = \/hD /UbT, where D is the diffusion constant. As 
a result, UCF has been observed only in ID (nanowires) 
and 2D systems (ultra thin-films and semiconductor de- 
vices) in samples with L < Lt ~ 1/im p3l . [2^ . Mag- 
netofingcrprint signals in large H have never been re- 
ported in mm-sized samples (in 2D or 3D). 

Hence, at T ^ 1 K, both the AB oscillations and 
UCF, which persist to intense -ff, are confined to /im- 
sized samples, whereas the AAS oscillations may be ob- 
served in mm-sized periodic arrays, provided H <'^200 
Oe. These comparisons show that the fingerprint sig- 
nal in Bi2Se3 is difficult to account for in terms of AB 
or AAS oscillations or UCF. Specifically, with D ^ 0.05 
m^/s, we find Lt ~ 1.5 /im at 1 K. The measured volume 
in G4 (2x10"'* cm^) exceeds the phase-coherent volume 
Lt by a factor of 6x10^. By classical averaging, the 
UCF should be strongly suppressed. Quantitatively, the 
scaling of the variance of G with L is derived |l13!] as 
Var\G{L)] - {e^ /hf{LT/L)^-'^, with d the dimension. 
This gives the rms amplitude of UCF as 0.01-0.05 /h, 
or 200-500 times weaker than depicted in Fig. [3l 

Further insight into the fingerprint is obtained from its 
autocorrelation function 

CiB, T) ^ , (1) 

which measures how a particular peak of 5G is correlated 
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FIG. 4: The magneto-fingerprint signal in tilted field. (Panel 
a) The MR curves in G4 at selected field-tilt angles 4° < 
61 < 85° at T = 2.2 K. Panel (b) shows the rms amplitude 
rms[5G] vs. tilt angle 9 (6 is defined in inset). The fit to 
rms[5G] = [a + h cos^ {e)\{e^ / h) yields a = 0.832 and h = 2.08 
(solid curve). Panel (c) plots the correlation function C{B) 
vs. H sX selected 9. 



with peaks at other field values {{■ ■ ■)b' means averaging 
with respect to i?', and B = fioH with /xq the perme- 
ability). As calculated [ll|, C(i3) for UCF should have a 
power-law decay ^ B'^~'^. 

In Fig. [2j;, we plot the correlation C{B, T) calculated 
from the results in Fig. [3^. Instead of a power-law de- 
cay, C{B) is oscillatory, displaying '-~^9 periods with an 
average period Bp '^1 T. If we interpret our period Bp 
as a characteristic length Lp ~ ^JW^JB^ ~ 640 A, the 
results suggest a characteristic area L^ in the network of 
conduction paths. 

The spin s of the carriers plays a significant role in the 
magnetofingerprint. We next describe the effect of tilting 
H at an angle 6* to c in the y-z plane. Figure [3^ shows 
plots of the conductance vs. H in Sample G4 at selected 
6*, with T fixed at 2.2 K. At each 0, the curve is comprised 
of the weak-field quantum anomaly superimposed on a 
parabolic (semi-classical) background. At 6' = 85°, the 
orbital term is negligible compared with the spin term 
because the flux inside the sample is reduced by a factor 
of 40. Hence the anomaly may be identified with the spin 
degrees alone. 

As in Fig. [31 we have extracted the conductance fluc- 
tuations at each 0. The gradual decrease of the ampli- 
tude rms[(5G'] vs. 6 fits well to [a -\- bcos^6]e'^/h (Fig. 
[4)3) . The first term a represents the spin term, while the 
second term (~ H^) is from orbital coupling. At maxi- 
mum tilt (85°), only the spin term a (comprising 29% of 
the rms) survives. This is rare example of a conductor 
in which the spin degrees are shown to generate a fin- 
gerprint signal. The correlation function C{B,6) of the 
fiuctuations (shown in Panel c) remains oscillatory with 
the same field period i3p ~ 1 T at all 9. 
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As displayed in Figs. [2^ and b, the weak-field anomaly 
in the MR has a \nH dependence that may be expressed 
as 

AG(iJ) = - ^ [Aorb + As,nn] In H, (2) 

where Aorb and Aspi„ (the orbital and spin terms, respec- 
tively) are both positive and comparable in magnitude. 
A positive Aorb implies field-suppresion of antilocaliza- 
tion in a 2D (two-dimensional) system (theory |25| pre- 
dicts Aorb — l/27r). In 2D systems, Coulomb interaction 
effects lead to a spin-Zeeman term of the form in Eq. 
[5] with A^n = F„/2-K, where the parameter F„ is of 
order 1 [I^. However, despite the suggestive InTJ de- 
pendence, we again encounter a large discrepancy in the 
size of the anomaly. The fit in Fig. yields Aorb = 122 
and Aspin = 78 for G4. Both are 500-700 x larger than 
the theoretical 2D values (see Table I). 

The most interesting question raised by these results 
concerns the nature of the states in the gap responsible 
for the magnetofingerprint signal in non- metallic Bi2Se3. 
Interpreting G in G4 as a 2D sheet conductance would 
give G/(e^//i) = kpC- ~8,000, which is far too large com- 
pared with kpl ^120 given by Sp measured in the metal- 
lic crystals. This implies that the in-gap states in G3-G8 



are bulk-like. On the other hand, the \ogH antilocal- 
ization anomaly identified in Eq. [2]) implies that they 
have 2D character (despite the large Aorb and Aspin)- 
We remark that there is evidence that ties the bulk in- 
gap states to the SS in 612863. In the APRES results, 
bulk states with 2D dispersion are seen to coexist with SS 
throughout the gap region (shown in yellow in Fig. la,b 
of Xia et al. They propose that the bulk states are 

confined by band-bending near the surface. If hybridiza- 
tion with the SS is strong, we expect that the bulk states 
will also display a large Rashba coupling that locks the 
spin transverse to k, thereby explaining the role played 
by the spins. If the 2D bulk gap states share the unusual 
properties of the SS, the long-range phase coherent na- 
ture of the fingerprint signal may be traced ultimately to 
the protected nature of the SS. The anomalies reported 
here provide a window into the unusual transport prop- 
erties of the SS. 
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